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A computational study of sliding blocks on inclined surfaces is presented. Assuming that 
the friction coefficient [I is a function of position, the probability -P(A) for the block 
to slide down over a length A is numerically calculated. Our results are consistent with 
recent experimental data suggesting a power-law distribution of events over a wide range 
of displacements when the chute angle is close to the critical one, and suggest that the 
variation of fj, along the surface is responsible for this. 

1. Introduction 

The dynamics of rigid bodies sliding on inclined planes under the action of 
gravity and possibly friction forces is a very old subject, present in most basic physics 
courses. However these systems are still being studied in order to better understand 
some of the intriguing properties of granular materials. Several recent papersB - ! 
present experimental and computational results that demonstrate the complexity of 
friction-related phenomena in specific cases. A recent study i of cylinders sliding on 
a rigid aluminum bar by Brito and Gomes shows that even this simple system may 
present unusual features. Measuring the number N(X) of slidings with length larger 
than A (see figure (1)) the authors find broad regions that can be very well fitted by 
a power law N(X)/N(Q) ~ \~ B , where the exponent B does not seem to depend on 
the material or the inclination 9 of the chute, when 9 is close to the critical angle 
9 C = arctan/i. Their data are in agreement with the Gutenberger-Richter lawu for 
the distribution of earthquakes and with numerical simulations by Chen, Bak and 
ObukhovB exhibiting self organized criticality. However, no theoretical explanation 
of this behavior was advanced. 

In this paper we present a simple model in order to explain the behavior observed 
in the above mentioned experiments. We consider the problem of a block that slides 
down, under the effect of gravity and friction forces, on an inclined surface at an 
angle 9 with the horizontal. We assume that friction is due to the existence of many 
uncorrelated contact points between the surfaces, and therefore becomes a rapidly 
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varying function //(x) of the block position x on the chute. At t = the block is set 
in motion with velocity vq. If 9 < 9 C , where 8 C = arctan/I is the "critical angle", 
the block will stop with probability one after some finite displacement A. We have 
in this case an 'avalanche' of size A. Our aim is to obtain the probability P(A) that 
the block stops at position A. 

The total energy variation after the slide is 

AE = mgAh - -mvl , (1) 

where Ah — hj — ho is the vertical displacement, the initial velocity and m the 
mass of the block. Since energy is dissipated by friction forces, we also have 

AE = / F^(x)dx=- mgcos(9){i(x)dx (2) 
Jo Jo 

where 9 is the inclination of the plane, x the block position along the plane, A(= Ax) 
the total displacement and fi(x) the local friction coefficient. Equations (Q) and 
(D) admit analytical treatment. The results of our theoretical approach will be 
published elsewhereO. 

2. The Computational Model 

In order to simulate the position-dependent behavior of \i we introduce a very 
simple computational model. This model is based on the following assumptions: 
We assume that friction occurs due to randomly scattered contact points between 
the two surfaces, and that these contact points are separated by a characteristic 
length a. The displacement of the block is discretized in small steps of length a, 
and we represent the "rugosity" of both surfaces by means of two binary strings of 
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Figure 2: An example of the binary representations of the block and chute. Only 
those regions contribute to friction for which both the block and the plane have a 
bit set to 1. 



Os and Is. Each bit can be thought of as representing the average properties of the 
surface over a length a. If for instance a certain region is more prominent than the 
average, the corresponding bit is set to 1, and to in the opposite case. Thus when 
the two surfaces are put in contact, only those regions will contribute to friction 
which both have a 1 on the corresponding location of their bit string. 

We model the disorder by assigning random strings of bits to both the plane and 
the block before each experiment. For computational convenience we set the block- 
string length to be 32 bits. Typical plane lengths are on the other hand 10 5 bits. 
The concentration of Is on the plane- and block-string are C p and Cb respectively. 

In Fig. (2) we show a schematic diagram of a block of length 14 bits with 
concentration Cb — 0.5 of ones, sliding over a plane with length L = 39 bits, with 
C p = 0.5. 

The local coefficient of friction ^i(x) is defined as 

/*(*) = (3) 

^ max 

where iV(a;) is the number of coincident Is and depends on both strings, N max is 
the block length in bits and 6 is a constant that can be associated to the strength 
of each individual contact. 

The block motion can be efficiently simulated in the following way: 

1. The block starts at x = with total energy E = \vhvq. 

2. We let the block slide over a distance a corresponding to one bit, and its 
energy variation is calculated as 
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AE = aF^(x) — mgasin(9) — mga cos(6) - sin(#) I (4) 

\ ^ max J 

3. If the total energy after this change turns out to be zero or negative, the block 
has stopped. Otherwise we set E — > E + AE and go to 2. 

The critical angle 9 C can be obtained by setting < AE >= in (||), and satisfies 

tan# c =~p = bCpCb (5) 




Figure 3: The accumulated distribution N(\) of slidings larger than A for several 
values of the inclination angle 9. Averages were done over 10 7 slidings. The critical 
angle 9 C is 7r/4. The straight line ^V(A) ~ A -1 / 2 shows the expected behavior for 
9 = 9 C . 



3. Results 

We measure the number N(\) of slidings with size larger than A. We fix C p — 
Cb = 0.5 and b = 4 so that 9 C — 7r/4. In Fig. (3) average results are shown for 10 7 
slidings on a plane of maximum length 10 5 bits and several values of 9 < 9 C . The 
straight line is shown for reference and corresponds to N(X)/N(0) = A~ 1//2 . Our 
numerical results indicate that, when 9 is close to 9 C , the distribution of avalanches 
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shows a power-law behavior with exponent 1/2 over a wide range of sizes. This 
exponent has the same value as obtained in experimental. As expected on simple 
grounds, below 9 C the distribution becomes exponentially decreasing for large sizes 
and therefore there is a finite average size. As already mentioned, this problem 
admits analytical treatment as well In Fig. (3) we show a preliminary theoretical 
result corresponding to (9 C — 9) = 10~ 4 . 

The average sliding size A was also measured, and the results are presented in 
Fig. (4). A power-law A ~ (9 C - BY" is obtained with v = 1.00 ± 0.02. This value 
is consistent with v = 1, which is obtained if /i does not depend on position. In this 

case A = 2 gsin°e -e) ~ $ c ~ f° r ^ dose to 9 C 




Figure 4: Mean sliding size as a function of inclination angle 9. We see that A ~ 
{9 C - 9Y V holds for 9 -> 9 C . 
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4. Conclusions 

A very simple model in which the coefficient of friction changes from point to 
point on the surface is able to reproduce a power-law behavior in the distribution of 
slidings of a block on an inclined chute, as recently observed in experiments. This 
holds for values of 8, the inclination angle, smaller than but very close to C . In this 
limit, the average sliding size diverges as (0 C — 

It is possible to do a theoretical study of this problem considering the random 
variation of the coefficient of friction. We show a preliminary result, and a complete 
study will be published elsewhere i. 
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